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I. INTRODUCTION 

In this paper, we explore the theoretical errors of 
simple expressions for oscillation probabilities found 
from approximate formulations for the purpose of 
analysis and prediction of experiment at current and 
proposed future neutrino facilities. Our assessment 
includes the interaction of neutrinos with electrons 
shown to be essential by Wolfenstein in 1978 til and 
later identified by Mikheyev and Smirnov [2| as a 
likely explanation of the deficit of solar neutrinos dis¬ 
covered experimentally by Davis a a. Subsequent 
to the work of Wolfenstein, the effect of matter on 
neutrino oscillations has been commonly explored 
using exact computer simulations. Computer simu¬ 
lations of neutrino oscillations in matter are found 
in Refs. [5] and references therein. Such numerical 
studies are preferred when density profiles with a 
spatial dependence are essential. 

Our assessment of expressions for the oscillation 
probability in the flavor-changing sector is 

made by comparing the approximate expressions to 
corresponding exact analytical expressions for the 
oscillation probability of 3 coupled Dirac neutri¬ 
nos. The expressions assessed include not only neu¬ 
trino oscillation probabilities found in Refs. [fl-Q but 
also neutrino oscillation probabilities from yet an¬ 
other approximate formulation we propose based on 
expressions for the oscillation probability found in 


our recently published paper [of evaluated with £- 
expanded eigenvalues. 

The exact analytical expressions for the oscillation 
probability used for the comparison are obtained 
from our recently-developed Hamiltonian formula¬ 
tion [g]. They are algebraic functions of the par- 
rameters of the the Standard Neutrino Model [lOj] 
(SNM), the neutrino eigenvalues, and include mat¬ 
ter effects. The expressions were obtained by using 
the Lagrange interpolation formula EH to exponen¬ 
tiate the neutrino Hamiltonian and have been used 
to verify our own exploratory studies [12-15J that 
used the expressions appearing in Refs. @-@|- 

In Sect. El we review the simple analytical results 
presently available and in Sect. ED we dis¬ 

cuss the assessment itself. Then, in Sect. eh we 
assess the oscillation probabilities found in Refs. @- 
@]. This assessment, which appears in Sect. EEEO 
confirms the finding of Ref. Q that the accuracy of 
the approximate analytic results of Refs. is lim¬ 
ited by the presence of branch points in the analytic 
structure of the eigenvalues of neutrinos propagating 
in matter. 

The analysis in Sect. IVIIII also suggests that the 
the simple expressions for neutrino oscillation prob¬ 
ability we find using our Hamiltonian formulation [9] 
are accurate to a few percent in all regions. It also 
identifies regions where the approximate analytical 
results are reliable by making numerical comparisons 
to exact results [Bi] . The approximate expressions we 
propose are conservatively estimated there to be ac¬ 
curate to a few percent within essentially all regions 
of interest. Approximate expressions of such accu¬ 
racy would, of course, obviate the need for exact 
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computer simulations under many circumstances. 

Combining the results of Sect. IVIII with numeri¬ 
cal results presented in Ref. !{|, we identify the re¬ 
gions where the approximate oscillation probabili¬ 
ties found from our proposed formulation are more 
accurate than those of Refs. (§-Q. In a future publi¬ 
cation [H], we show how to construct such approx¬ 
imate oscillation probabilities in practice guided by 
the assessment made here. 


II. NEUTRINO DYNAMICS 

The dynamics of the three known neutrinos and 
their corresponding anti-neutrinos in matter is de¬ 
termined by the time-dependent Schroedinger equa¬ 
tion, 

> = H v \v{t) > , (1) 

where the neutrino Hamiltonian H u , 

H v = H 0v + H x , (2) 

consists of a piece Hq v describing neutrinos in the 
vacuum and a piece H i describing their interaction 
with matter. 

The solutions of Eq. CD may be expressed in terms 
of stationary-state solutions of the eigenvalue (EV) 
equation 


depends on time only through the time difference 
t' — t. Then, written in terms of the stationary state 
solutions | v m i > of Eq. 0, 

= ^2 I "mi > e ~ lEi(t ~ l) < v mi \ . (5) 


We assume here that that neutrinos and anti¬ 
neutrinos represented by |rW > and | v® ni >, respec¬ 
tively, are the structureless elementary Dirac fields 
of the the Standard Neutrino Model jlc|. For this 
reason the theory is invariant under CPT, so the 
mass of an anti-neutrino in the vacuum is the same 
as that for its corresponding neutrino. 

The case of main interest for many situations is 
the ultra-relativistic limit, \p\ » m 2 (we take the 
speed of light c = 1). For ultra-relativistic neutrinos 
in the laboratory frame, the energy of a neutrino in 
the vacuum becomes 




( 6 ) 


where nii is its mass the vacuum. Similarly, Ei ap¬ 
pearing in Eq. ([3]) may be written 


M 2 


m 


where Mi is its mass in the medium. Thus, in this 
limit and in dimensionless variables, 


Hv\Vmi ^ > > ( 3 ) 

where the label “to” indicates neutrino mass eigen¬ 
states, as distinguished from their flavor states some¬ 
times denoted the label “/”. Neutrinos are produced 
and detected in states of good flavor. 

In operator form, the dynamics of neutrinos may 
be expressed in terms of the time-evolution operator 
S(t',t), which describes completely the evolution of 
states from time t to t' and also satisfies the time- 
dependent Schroedinger equation. 

We will examine neutrinos propagating in a uni¬ 
form medium for interactions constant not only in 
space but also time. Because the Hamiltonian is 
then translationally invariant, attention may be re¬ 
stricted to states, both in the vacuum and in matter, 
characterized by momentum p and therefore having 
the overall r-dependence e lp ' r . In this case, expres¬ 
sions may be simplified by suppressing the overall 
plane wave, a convention we adopt. 

For time-independent interactions, ), 

S{t',t) = e- iH ^ t '- t) , ( 4 ) 


and 


with 




M 2 — m\ 

9 9 " 

771 g — 7711 


/ ° o ° \ 

Hqv —> I 0 a 0 ] 

V° 0 


2 2 

777-2 — m l 

2 2 ’ 

TTIg — 771 ^ 


( 8 ) 


(9) 


( 10 ) 


In this limit, the distance L from the source to the 
detector corresponding to S(t',t) in Eq. Q is 

L = t'-t. (11) 

The time-evolution operator, Eq. 0, expressed in 
dimensionless variables is then, 

S(L) = 

_ g2iE 1 Aig— 2i£[„Al 


( 12 ) 




3 


where H u is the full neutrino Hamiltonian expressed 
in dimensionless variables, and where 


L{m\ - m\) 
L ~ 4 E 


(13) 


III. THE NEUTRINO INTERACTION AND 
THE STANDARD NEUTRINO MODEL 

Flavor and mss states are related by the neutrino 
analog of the familiar CKM unitary matrix U, 

Vf = Uv m . ( 14 ) 

The matrix U is often parametrized in terms of three 
mixing angles (9 12, 9 13 , 023) and a phase 5 cp charac¬ 
terizing CP violation, taking the form, 

/ C 12 C 13 S 12 C 13 si 3 e- <4 <* \ 

U = I U 21 U 22 S 23 C 13 I ; (15) 

\ U31 U32 C23C13 J 

where, 

U21 = — S12C23 — Cl 2 S 23 Sl 3 e i,5cp 

U22 = C12C23 — Sl 2 S 23 Sl 3 e l,5cp 

U31 = S12S23 — Ci 2 C 23 Sl 3 e* 5cp 

U32 = —C12S23 — Si 2 C 23 Sl 3 e i<5cp , ( 16 ) 

and where the index i = 1 in Uij corresponds to 
the electron (e) neutrino, i = 2 to the muon (p) 
neutrino, and i — 3 to the tau (r) neutrino. The 
standard abbreviations, S 12 = sin$ 12 , C 12 = cos $ 12 , 
etc have been used. 

The perturbing Hamiltonian H 1 is determined 
by the interaction between electron neutrino flavor 
states and the electrons of the medium. Expressed 
as matrix, 

(V 0 0 \ 

H 1 = U- 1 \ 0 0 OP, ( 17 ) 

\ 0 0 0 / 

with V = ±\/ 2 Gpn e and n e the electron number 
density in matter. 

For electrically neutral matter consisting of pro¬ 
tons, neutrons, and electrons, the electron density 
n e is the same as the proton density n p , 

Tt e = Tip 

= RA , ( 18 ) 

where A = N + Z is the average total nucleon num¬ 
ber density and R = Z/A is the average proton- 


nucleon ratio. In the earth’s mantle, the domi¬ 
nant constituents of matter are the light elements 
so R ~ 1/2; in the surface of a neutron star R « 1. 
Matrix elements of H\ are thus 

< M{k)\Hi\M{k') > = U* lk VU lk , . (19) 

Using the well-known expression for V, we find 
the corresponding the (dimensionless) interaction 
strength A of neutrinos and anti-neutrinos with mat¬ 
ter to be, 

A = ±6.50 10 ~ 2 R U[GeV]p[gm/cm 3 ] , ( 20 ) 

with E[GeV] being the neutrino beam energy E 
(in GeV) and p[gm/cm 3 ] the average total density 
(in gm/cm 3 ) of matter through which the neutrino 
beam passes on its way to the detector. For our 
experiments close to the earth’s surface, the appro¬ 
priate density is mean density of the earth’s mantle, 

p [gm/cm 3 ] = p 0 

« 3 . ( 21 ) 

We adopt the SNM 0 , given next, to complete 
our description of the neutrino Hamiltonian. Most 
of the parameters of the SNM are consistent with 
global fits to neutrino oscillation data with relatively 
good precision Emu. These include the neutrino 
mass differences, 

2 2 _ c 2 

m 2 — m 1 = om 2 i 

= 7.6 x 1CT 5 eV 2 (22) 


and 


2 2 _ c- 2 

m 3 — m 1 = om 31 

= 2.4 x 1(T 3 eV 2 , (23) 


which corresponds to 

= Sml 1 
^ _ ^31 

= 3.17 x 1CT 2 . (24) 

The mixing angles 9 are also determined from ex¬ 
periment. In the SNM, the value of $ 23 > 

023 = tt /4 , ( 25 ) 

is the best-fit value from Ref. [Tj, and $ 12 , 

0 12 = tt/5.4 , (26) 

is consistent with the recent analysis of Ref. Il3 |. 
The mixing angle $13 is known to be small ($13 < 




4 


0.18 at the 95% confidence level) but until recently 
its precise value is uncertain. Results from the Daya 
Bay project 51 have measured its value quite acu- 
rartely, sin $13 « 0.15, which we adopt to determine 
our value for $ 13 , 

0 13 = 0.151 . (27) 

This fixes R p = sin 2 $ 13 / 0 : sa 0.711. The CP vio¬ 
lating phase is not known at all and will one of the 
major interests at future neutrino facilities. 

Using the value of dm\ x from the SNM, A^, de¬ 
fined in Eq. (1131) becomes, in the high-energy limit, 

“ 3 05 * - < 28 > 

Here L [Km] is the baseline and E[GeV] is the neu¬ 
trino beam energy 

IV. EXACT ANALYTIC EXPRESSIONS 

FOR V(yi —>■ Vf) 

Our Hamiltonian formulation formulation of neu¬ 
trino oscillations [9] leads to exact, closed-form, an¬ 
alytical expressions neutrino oscillation probability 
expressed as the sum of four partial oscillation prob¬ 
abilities, 

V(v a —> V b ) = 6(a , b) + P s “n ,5 + Po b + PcosS 

+ P c a o S *6, ( 29 ) 

representing the dependence of the oscillation prob¬ 
ability on the CP violating phase 5 cp . The partial 
oscillation probabilities are given in Appendix [A] 
where they are expressed as explicit functions of the 
neutrino medium modified neutrino eigenvalues Eg 
and a set of coefficients w°' b f: . 

Equations determining the exact eigenvalues, and 
various approximations to them, appear in Ap¬ 
pendix IA II Finally, in Appendix E explicit alge¬ 
braic expressions for the coefficients for all tran¬ 
sitions v a tA are given in terms of the parameters 
of the SNM. 


V. PERTURBATIVE EXPANSIONS OF 

V(vi -A v f ) 

The fact that a and sin $13 are naturally small 
in the SNM commonly motivates approximation 
schemes |2C} based on first-order Taylor se¬ 


ries expansions in one of these small parameters, 
(where ^(stands for a or sin 2 $ 13 ). For example, in 
Refs. |3, @1 the oscillation probability is expanded 
in sin 2 $ 13 . Reference Q makes use of an expansion 
in the small parameter a. In both publications, the 
only transitions considered are those in the (e,/z) 
sector. 

Although the expansions have the advantage of 
simplifying the theory, their use comes at a price: 
neither expansion gives an accurate representations 
for all regions of the interaction strength A 1 in¬ 
cluding some regions of critical importance. The 
origin of this inaccuracy is the presence of branch 
points in the neutrino eigenvalues 0. We refer to 
these truncated ^-expanded oscillation probabilities 
of Refs, m as “full” results. In addition to his full 
result, Freund also introduces an ad hoc “patched” 
result to repair spuriosity of different origin. 

In this section, we review the results for the os¬ 
cillation probability in the (e, /z) sector obtained in 
Refs. @-@, [ 20 }. In later sections, we then com¬ 
pare their oscillation probabilities numerically to the 
corresponding exact expressions obtained from our 
Hamiltonian formulation evaluated with both the 
exact neutrino eigenvalues and with the ^-expanded 
eigenvalues. 


A. Full Result of Freund 

In his work, Freund gives @ two versions of the 
oscillation probabilities based on this a-expansion, 
a full version and a patched version. As in our own 
work, the oscillation probabilities are expressed as a 
functional of the flavor mixing angles and the CP- 
violating phase that define the neutrino CKM mix¬ 
ing matrix, the neutrino mass differences squared, 
and the strength of the neutrino interaction with 
matter. Consequently, Freund makes simplifica¬ 
tions [fi] in two stages. The first stage is to expand 
the eigenvalues, and the second stage is to use these 
results in an expansion of the oscillation probability. 

Freund’s a-expanded eigenvalues, given in 
Eqs. (18,19) of Ref. @ are identical to those of 
Eq. (??) in the Appendix. Numerical comparison to 

~ a. ±cx 

the exact result confirms that E 1 and E 2 are poor 
representations of the corresponding exact results 
in the vicinity of the branch point at A = 0. 

Freund’s full a-expanded partial oscillation prob¬ 
abilities in the v e —> v p sector are, 
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P^(A l ,A) = sin 2 0 23 Sm 2 f 13 sin 2 A L C a 

ct 

r)GU / a a x ^ ^ 2013 . 2 d • 2 a/i *2/i a 

P 1 ^(A L ,A) = -a - -p- -sm 0 i2 sm 20 i3 sm 9 23 A L 


Cl 

x sin2(7 a AL + — ^4 +^cos.. 0 i 3 ) ^ g j n 2 2 g 13 s j n 2 g 23 s j n 2 A^Ca 

6 »i 3 sin 2012 sin 26>i3 sin 26>23 . a a 

- — - sm Cv»Ar 

\r< a 


a . cos 0 i 3 f 


P e P cfAi, A) = —sin< 5 C r,-—- 

sm5V ’ 2 cp AC a cos^ 0 i3 

x [cosCqAl — cos (1 + A)Al\ 

a cos 0i3 sin 2012 sin 2013 sin 2023 A A 

cos 8y^-L,si) = — cos o cp -—-—-SinAiCi 


P cos6( A L,A) = ^COsS cp - jp. -—- 

^ AC a COS- 013 


- y Ct 1 

iC a A L ] 


x [sin (1 + A)Al =F sin C a 

n e WA A^ Tl + Ca ± A COS 2013 n . nn . n/1 

P 2 (Al, A) = a -- —7 -cos 0 i 3 sm 2012 sm 2013 

2C%A cos 2 0 i 3 

x sin 2023 sin 2 C a A^ 

rpefj, / a tx 2 2(7,, COS 2 023 sin 2 2012 . 2 1 , i ^ A x a 

P^{A L ,A)=a -j- - ' ~ - — - sin -(1 + A T C a )A L 

A 2 cos 2 Oi 3 (^fA + C a ± cos 20 13 ) ^ 


(30) 


The lower sign applies above the atmospheric reso¬ 
nance, and the upper sign applies below it. The term 
P 3 , of order a 2 , is a correction term that improves 
the overall accuracy. To obtain these results, Fre¬ 
und linearized the oscillating terms over a, the solar 
mass-squared difference. His results thus apply only 
for baselines where 

aA L < 1 . (31) 

In addition to observing that his expressions for 
the a-expanded eigenvalues are not valid for |x4| < a , 
Freund [6j notes a shortcoming of the technique he 
uses when when there are resonances. The advan¬ 
tage of this technique is the possibility of extract¬ 
ing the mixing angles, as modified by the medium, 
from the unitary transformation that diagonalizes 
the Hamiltonian. 

However, in order to construct this transforma¬ 
tion, it is necessary to keep track of the ordering of 
the eigenvectors across resonances that characterize 
the transition of interest. In the presence of reso¬ 
nances, the eigenvectors must be exchanged at their 
position to preserve the proper ordering. Unavoid¬ 
ably, doing so introduces discontinuities. This is the 
source of the flip of sign at at the location of the 


atmospheric resonance, A = cos20i3, in Eq. ®. 
This discontinuity is the price one pays to use this 
method. We note in passing that these discontinu¬ 
ities are avoided in our Hamiltonian approach @ by 
using the Lagrange interpolation formula. 

Freund examines the accuracy of the result given 
in Eq. (1301) by comparing to exact numerical sim¬ 
ulations for a specific baseline and range of ener¬ 
gies encompassing the atmospheric resonance. From 
these results, given in his Fig. 1, he concludes that 
the loss of precision of his full a-expanded is less 
serious for small values of 0 i 3 , such as the value 
013 ~ 0.1, which is the CHOOZ bound. Unfortu¬ 
nately, the recent result Daya Bay JT^| gives a larger 
result, sin 0i3 ss 0.15, which makes the use of Fre¬ 
und’s full result somewhat awkward. 


B. Patched Result of Freund 

Freund identifies the source of the inaccuracies to 
be the sub-leading terms in 0 i 3 within his representa¬ 
tion. Based on his understanding, Freund proposes 
a “patched” result that he argues would be more 
appropriate than his full result when it fails. This 
leads to his “patched” expressions, 
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Po M (AL,i) 

P s Zs(^l,A) 

p£ s (a l ,A) 


■ 2 a s i n2 2013 . 2 / a \ A 
sin 0 23—-—sin (A-l)Ax 


t sin 5. 


cp 


; COS <5 C 


(A — l ) 2 
cos 0 i 3 sin 2012 sin 20 i 3 sin 20 2 3 

A(A-1) 

cos 0 i 3 sin 2012 sin 2013 sin 20 23 

1(1^7) 


sin Al sin AA x sin (.A — 1 ) Ax, 
cos Ax sin AAl sin (A — 1 )Ax 


P^(A L ,A) = a 


2 COS 2 023 sin 2 2012 . 2 


A 2 


sin 2 AA l 


(32) 


which appear in Eq. (38a-38d) of Ref. Freund 
notes that the 0 i 3 dependence is improved by re¬ 
taining P^(A l ,A). He neglects P^' i (A L ,A) and 
P| M (Ax,A) because he finds them much smaller 
than PZ s (Ax, A) and P^ s 5 (Ax, i). 


C. Full Result of AHLO 


The AHLO neutrino oscillation probability, 

V(v a -)■ v b ) = P v 

= |S ab | 2 , (33) 

is determined from an expansion of time-evolution 
operator S ab in sin 2 0i 3 . The approach was origi¬ 
nally proposed in Refs [7j, , and details of its im¬ 

plementation are given in Appendix A of In Ref. Q. 
Their final result is a first-order perturbation theory 
expansion of S ab in sin 2 0i 3 , denoted there S'(t,to) 
and is equivalent to a first-order Taylor expansion of 
the entire V(yi —> v/) about = sin 2 0i 3 = 0, 

P(v a —> v b ) re P a6 || i=sin 2 g 13=0 

d p ab 

+ sin 2 0i 3 —|{ i= o • (34) 

i 

Because this expansion involves the total derivative, 
dP ab /dP, the expanded eigenvalues are not used ex¬ 
plicitly El- 

The expression for S'(t,to) is found from the so¬ 
lution of the time-dependent Schroedinger equation 
given in Eq. (A.9) of Ref. Q, 

i—Sq^Ao) = HqSq^, to) , (35) 

with the Hamiltonian H' 0 , expressed in the neutrino 


mass basis given in Eq. (A. 6 ), 


—a cos 2012 + A a sin 2012 
a sin 2012 acos 20 12 — A 


,(36) 


where 


A = 


2 E 


(37) 


The solution of Eq. (1551) . Eq. (22) of Ref. 1, 


So(t> to) 


a (Mo) P(t,to) \ 
-P*(t,to) a* (Mo) ) ’ 


(38) 


is expressed in terms of a(t,to) and /3(Mo)- The 
quantities a(t,to) and /3 (Mo)j which may be found 
in closed form, determine S ab to leading order in 
sin 0 i 3 , 


s '(t, t 0 ) — ( s 21 • (39) 

All higher-order corrections are contained in the 
transformation that rotates S' 0 back to the original 
flavor basis. 


Defining the baseline L in terms of the time t — to 
for the neutrinos to travel between their produc¬ 
tion target and detector in a uniform medium, as 
in Eq. CD. all quantities become functions of L. In 
particular, £i 2 (Mo) -A S\ 2 (L), where 

£12 (T) = cos 0 23 /3(T) - i sin 0 23 A(L) . (40) 

The quantity A(L) = aA a (L) + bA b (L) is given as 
the first expression in Eq. (A. 14) of Ref. [7j, with 
A a (L) defined in Eq. (A.16) and A b (L) in Eq. (A.17). 
These quantities may be calculated analytically as 
integrals over L involving f(L), 

f(L) = e" iAi , 


(41) 









7 


a(L), and /3 (L). We have solved these equations 
for S'(t 7 to) exactly, and all numerical AHLO results 
are obtained from these expressions. These result¬ 
ing expressions are, however, is very complicated 
and too awkward to use in practice. We have used 
them because our only interest here is to benchmark 
Refs [3, H[ . Simplified expressions for the matrix ele¬ 
ment S 12 , which describes —> v e transitions, were 
used in Refs. fl2l - fl5l |. Because these expressions en¬ 
tail a yet further set of approximations, for the pur¬ 
pose of this paper it serves no purpose to examine 
these simplified expressions. 

Below, we numerically compare the oscillation 
probability obtained from S'(t,to) to the exact one 
obtained from our Hamiltonian formulation. The 
comparison is made for the specific case of e —> p 
transitions. 


VI. ASSESSING THE ^-EXPANDED 

V(ui ->■ u f ) 

In Sect. IVII1 we compare the a-expanded oscilla¬ 
tion probabilities, the sin 2 0 13 -expanded oscillation 
probabilities, and the exact oscillation probabilities. 
From these results, we are able to identify the re¬ 
gions where both ^-expanded oscillation probabil¬ 
ities might be improved by using our Hamiltonian 
formulation. 


A. Regions of A 

We assess the expanded oscillation probabilities 
over specific regions of A and A L . For this purpose, 
the interval A > —0.8, broken down into five sub- 
regions spanning the interval Ai < A < Af. 

The first region is the solar resonance region con¬ 
taining the solar resonance at A = a. This region 
extends from A; = 0 to Af = 0 . 1 . Motivated by 
the representation of the sin 2 0 i 3 -expanded eigenval¬ 
ues given in Appendix IA 11 we split the solar reso¬ 
nance region into the deep solar resonance region, 
0 < A < a and the far solar resonance region, 

a < A < 0 . 1 . 

The second region, the transition region, extends 
from Ai = 0.1 to Af = A 2 , where A 2 ~ 0.538 
is the value of A at which the periods of the two 
most slowly varying Bessel functions, jo(A[l]) and 
jo (A [3]), are equal. The third region includes the 
atmospheric resonance at A = Ao ~ 1 and extends 
from Ai = A 2 , Af = 1.2, The asymptotic region in¬ 


cludes all values of A > 1 . 2 . Finally, we consider the 
region extending from Ai = —0.8 to Af = +0.8. 

Within each region considered, is chosen so 
that 'P(r'e —> I'iA displays the greatest sensitivity to 
approximations. This is assured by requiring that 
each of the three Bessel functions constituting the 
oscillation probability are of comparable size and in¬ 
terfere both constructively and destructively. For 
maximal interference, each of these Bessel functions 
jo(A^) should undergo at least one period of oscilla¬ 
tion when the neutrino is detected at a distance A l 
from the source. This is fully discussed in Ref. [&]. 

The maximal interference within the solar reso¬ 
nance region was determined to occur for A^ = 60; 
the transition region for Al = 17; the atmospheric 
resonance region, for Al = 4; the asymptotic re¬ 
gion for A/, = 4. For the region extending from 

Ai = - 0.8 to A f = + 0 . 8 , a l = 35 . 

The Ai dependence is examined for two values 
of A, one in the solar resonance region and one in 
the atmospheric resonance region. In all cases, the 
oscillation probability is examined for S cp = 7r/4, 
with the other parameters taken from the SNM. 


B. Extrapolations 

The results shown in this paper may be extrapo¬ 
lated to a variety of values baselines, medium prop¬ 
erties, and neutrino energies within the regions un¬ 
der using Eqs. (1131201) . 

r n 3 I Al A^ 

L Km = 5.04 x 10 3 - 1 1 „ 

l?p[gm/cm 3 ] 

= 15.3-^^ , (42) 

i?p[gm/cm 3 ] 

for a medium described by p and R. Using Eq. (H^l) . 
the oscillation probability may be extrapolated to a 
variety of values of neutrino energy, base lines, and 
medium properties once it is known for a specific set 
of values (A l, A) 



VII. COMPARISON OF THE FULL AND 
PATCHED V{u e -> rv) OF FREUND TO THE 
FULL V{v e u,,) OF AHLO 


In this section, we assess simplified expressions for 
V(y e —» i+) found in Refs pial- These expressions 
are the only published attempts made to simplify 
the oscillation probability. 








We compare the simplified P(i/ e —► uA to the cor¬ 
responding exact oscillation probability [9J evaluated 
with the exact eigenvalues to determine the relative 
and absolute accuracy of the approximate oscilla¬ 
tion probabilities commonly used for the analysis 
and prediction of neutrino oscillation phenomena. 
From these results, one easily identifies which of the 
£i-expanded oscillation probabilities best describes 
each region. 

Freund’s patched is shown as a short-dashed 
curve, his full result as a medium-dashed curve, the 
AHLO oscillation probability as a dot-dashed curve, 
and the oscillation probability of our Hamiltonian 
formulation Q evaluated with the exact eigenvalues 
as a solid curve. 


A. A l Dependence of the ^-Expanded 

V(u e ->■ Vp) 

In this section, we compare the oscillation proba¬ 
bilities over intervals of for several values of A, 
one near the solar resonance and the other near the 
atmospheric resonance. 

1. Near the Solar Resonance 

We first look at the A^-dependence of 
P e_ ^(A^, A) in the solar resonance region for 
A = 0.0102. Figure [T] shows that Freund’s full 
a-expanded oscillation probability and his patched 
oscillation probability are both in reasonable agree¬ 
ment with the exact result for Al < 20, at which 
point both of Freund’s oscillation probabilities begin 
to depart from the exact oscillation probability. 
They increasingly diverge from the exact result as 
Al increases, which reflects the breakdown of the 
a-expansion in the solar resonance region. On the 
other hand, the AHLO result agrees with the exact 
result comparatively well over the entire interval of 
Al shown. 

2. Near the Atmospheric Resonance 

We next look at the A^-dependence of 
P e ^ rfJ -(A L , A) in the atmospheric resonance region 
for A = 0.8, a value of A below the atmospheric 
resonance. Figure [2] shows that Freund’s full 
a-expanded oscillation probability is in reasonable 
agreement with the exact result over the entire 
interval of A^ shown. On the other hand, Freund’s 





FIG. 1 . P efl (Al, A) in the solar-resonance region (A = 
0.0102) over the interval 5 < Al < 55 for neutrinos 
in matter. Parameters are taken from the SNM. Exact 
result (solid curve). Patched a-expanded result of Fre¬ 
und [ff (short-dashed curve). Full a-expanded result of 
Freund [ 6 | (medium-dashed curve). Full sin 2 $13 result 
of AHLO [?] (dot-dashed curve) 


pep 



FIG. 2 . Freund P bij, (Al,A) in the atmospheric reso¬ 
nance region (A = 0 . 8 ) over the interval 0 < Al <35 
for neutrinos in matter. Parameters are taken from the 
SNM. Exact result (solid curve). Patched a-expanded 
result of Freund [t| (short-dashed curve). Full a- 
expanded result of Freund 161 (medium-dashed curve). 
Full sin 2 #13 result of AHLO [7'] (dot-dashed curve) 


patched a-expanded oscillation probability and the 
AHLO result are quite similar but are in substantial 
disagreement with the exact result over the same 
interval of A^. 


B. The A Dependence of the ^-Expanded 

V(ui v f ) 

Below, we examine the A dependence over the five 
regions defined earlier. 
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FIG. 3. P e,J (A l, A) over the deep solar resonance re¬ 
gion, 0 < A < a, with Ai, = 60 for neutrinos in matter. 
Parameters are taken from the SNM. Exact result (solid 
curve). Patched a-expanded result of Freund [g] (short- 
dashed curve). Full a-expanded result of Freund 6 ] 
(medium-dashed curve). Full sin 2 #13 result of AHLO 31 
(dot-dashed curve) 


FIG. 4. P £IJ '(Al, A) over the far solar resonance region, 
a < A < 0.1 with Al = 60 for neutrinos in matter. 
Parameters are taken from the SNM. Exact result (solid 
curve). Patched a-expanded result of Freund [g| (short- 
dashed curve). Full a-expanded result of Freund 6] 
(medium-dashed curve). Full sin 2 813 result of AHLO 3 
(dot-dashed curve) 


1. Solar Resonance Region: 0 < A < 0.1 

The deep solar resonance region, which extends 
over the interval 0 < A < a, is shown in Fig. [3] 
The far solar resonance region, which extends over 
the interval a < A < 0.1, is shown in Fig. [4] The 
results of Freund Q and AHLO 3] are compared 
to the exact result within these two regions. The a- 
expansion is known to fail within the solar resonance 
region. Q because one of the eigenvalue branch point 
lies at a = 0. 

Although Freund asserts that his results fails only 
below the solar resonance, A < a, it is clear from 
Fig. 0 that the a-expanded results are poor approx¬ 
imation out to A ft 0.1. The failure becomes in¬ 
creasingly severe as a —> 0. This is evident for both 
Freund’s patched and full results. 

Clearly, the full sin 2 613 -expanded result of AHLO 
does much better within the solar resonance region. 
However, Fig. [3] shows that the AHLO result does 
not approach the exact result in the limit a —1 0. 


2 . 0.1 < A < A 2 



FIG. 5. P em (Al,A) for Al = 17 over the over the 
transition region, 0.1 < A < A 2 , for neutrinos in mat¬ 
ter. Parameters are taken from the SNM. Exact result 
(solid curve). Patched a-expanded result of Freund [§] 
(short-dashed curve). Full a-expanded result of Fre¬ 
und )6| (medium-dashed curve). Full sin 2 813 result of 
AHLO [3 (dot-dashed curve) 

transition region. 

Comparing the dot-dashed and solid curves, it is 
clear that that Freund’s patched and the AHOL re¬ 
sult provide rather similar descriptions of the oscil¬ 
lation probability throughout the transition region. 


In Fig. O we compare the oscillation probabilities 
within the transition region. Comparing the short- 
dashed and medium-dashed curves to the the solid 
curve, it is clear that that Freund’s patched and full 
results agree well with the exact result for 0.1 < 
A < 0.3, but for larger values of A his patched result 
increasingly departs from the exact result, whereas 
his full result tracks the exact result throughout the 


3. A 2 < A < 1.2 

In Fig. [Gl we compare the oscillation probabili¬ 
ties over the atmospheric resonance region. Com¬ 
paring the dot-dashed and medium-dashed curves, 
it is clear that that the departure of Freund’s full 
result from his patched result continues to increase 
across the atmospheric resonance region. Compar- 
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FIG. 6 . P eM (Az,, A) for A l = 4 over the atmospheric 
resonance region, A2 < A < 1 . 2 , for neutrinos in mat¬ 
ter. Parameters are taken from the SNM. Exact result 
(solid curve). Patched a-expanded result of Freund [gj 
(short-dashed curve). Full a-expanded result of Fre¬ 
und i (medium-dashed curve). Full sin 2 #13 result of 
AHLO 3 (dot-dashed curve) 


FIG. 7. P e,j (Al, A) vs A for Al = 4 within the asymp¬ 
totic region, 1.2 < A < 2.5 for neutrinos in matter. Pa¬ 
rameters are taken from the SNM. Exact result (solid 
curve). Patched a-expanded result of Freund [g| (short- 
dashed curve). Full a-expanded result of Freund 6 ] 
(medium-dashed curve). Full sin 2 613 result of AHLO 7| 
(dot-dashed curve) 


ing the dot-dashed and solid curves, it is clear that 
that Freund’s patched and the AHOL result con¬ 
tinue to provide rather similar descriptions of the 
oscillation probability throughout the atmospheric 
resonance region. 

Freund’s full result provides a reasonable descrip¬ 
tion of the exact result below the atmospheric reso¬ 
nance, at which point his result is discontinuous and 
lies above the exact result for larger values of A. 


4 . 1.2 < A < 2.5 

In Fig. [3 we compare the oscillation probabilities 
over a portion of the asymptotic region. Comparing 
the dot-dashed and short-dashed curves, we see once 
again the departure of Freund’s patched result from 
his patched result. However, this departure appears 
to decreases above the atmospheric resonance region 
and approaches the exact result for the larger values 
of A. 


5. -0.8 < A < 0.8 

Because the branch point of the a expansion oc¬ 
curs for A = 0, it will adversely impact the total 
oscillation probability for both neutrinos and anti¬ 
neutrinos This is confirmed in Fig. [5] 



FIG. 8 . P e>1 (Al, A) over the interval —0.8 < A < 0.8 
with Al = 35 for neutrinos in matter. Parameters are 
taken from the SNM. Exact result (solid curve). Patched 
a-expanded result of Freund [ 6 ] (short-dashed curve). 
Full a-expanded result of Freund [gj (medium-dashed 
curve). Full sin 2 $13 result of AHLO 0] (dot-dashed 
curve) 


C. Discussion 

In this section, we calibrated the accuracy of Fre¬ 
und’s 0 full and patched exact oscillation probabil¬ 
ities and the full AHLO [ij oscillation probability. 
We were able to identify which of the approximate 
solutions best described the exact result within and 
above the solar resonance region. Few studies of this 
nature have ever been made. 

Freund’s patched result was found to be signif¬ 
icantly less accurate than the full result in all re¬ 
gions. Freund’s full result was calculated and com¬ 
pared to the exact result in all regions, including the 
solar resonance region. The failure of Freund’s full 
result within the solar resonance region is evident 
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from Figs. 13141 we also confirmed that the sin 2 $13 
expansion leads to inaccuracy in the vicinity of the 
atmospheric resonance. Both inaccuracies grow as 
A approaches the branch points. 

In the vicinity of the atmospheric resonance, A « 
cos 2013, we found that neither the AHLO result nor 
the full result of Freund provides a very good de¬ 
scription of the oscillation probability. This is ex¬ 
pected for the sin 2 $13 expansion because the corre¬ 
sponding branch point lies close to the atmospheric 
resonance. Freund’s full result is discontinuous at 
the position of the atmospheric resonance. This dis¬ 
continuity is quite noticeable when the Daya Bay 
value of 0 i 3 is adopted. 

It was clear from our calculations that our Hamil¬ 
tonian formulation needs no patch to avoid any dis¬ 
continuity such as that present in Freund’s full re¬ 
sult. 

Freund, in his paper [g], does actually compare 
his full and patched results to an exact calculation 
of the oscillation probability [his Fig. (1)]. His re¬ 
sults, when viewed casually, may appear inconsistent 
with those given here, especially with regard to the 
theoretical error of his patched result. 

However, one should note that Freund’s compar¬ 
isons are presented on a logarithmic plot of multi¬ 
ple decades, whereas ours are given on a linear plot. 
On a logarithmic scale, factors of two are often get 
lost in the noise. Additionally, Freund considers a 
broad range of model parameters, which makes it 
difficult to pinpoint the portion of his presentation 
that reflects the SNM, which is of course the model 
underlying our analysis. 

We also found that the sin 2 0 i 3 expansion 3, [g] 
represents of the exact oscillation probability rela¬ 
tively accurately within the solar resonance region, 
but its accuracy in other regions is comparable to 
Freund’s patched result. We also found that in 
the region between the solar and atmospheric res¬ 
onances, 0.1 < A < A 2 , the expansion in a is a good 
representation of the exact result, whereas the ex¬ 
pansion in sin 2 0 i 3 rapidly deteriorates with increas¬ 
ing A. 

In the vicinity of the atmospheric resonance, A ss 
cos 2013, we found that neither the AHLO result nor 
the full Freund result provides a very good descrip¬ 
tion of the oscillation probability. This is expected 
for the sin 2 0 i 3 expansion because the correspond¬ 
ing branch point lies close to the atmospheric reso¬ 
nance. A careful reading of Ref. jg] reveals that some 
inaccuracy of the a expansion would have been an¬ 
ticipated there in light of the Daya Bay result pjj, 
but its extent may be a surprising disappointment 


to some. The asymptotic region, A > 2.0 is shown 
in Fig. □ 

Even though expanded eigenvalues are not used 
explicitly for the sin 2 0 i 3 -expanded oscillation prob¬ 
ability, Eq. (l34l) is still a poor representation of 
the exact result in the vicinity of the branch point 
at A = Aq, where the eigenvalue expansion fails. 
To separate the theoretical errors into those aris¬ 
ing from the expansion of the oscillation probability 
from those arising from the expansion of the eigen¬ 
values we have relied on the numerical techniques. 


VIII. APPROXIMATE FORMULATION OF 
P{ui ->■ u f ) ACCURATE TO A FEW PERCENT 

Because subtle but important trends are more eas¬ 
ily identified by examining analytic expressions than 
by scanning through numerical lists, simpler and 
most accurate analytic expressions would be signif¬ 
icant considering that higher quality data are ex¬ 
pected to become available at future neutrino facili¬ 
ties. In our original publication Q, we compare ex¬ 
act oscillation probabilities to those found from our 
Hamiltonian formulation evaluated with ^-expanded 
eigenvalues. 

A. Above the Solar Resonance Region, A > 0.1 

In our original publication [ 9 }, we find that the er¬ 
rors of 0 (a 2 ) expected in regions sufficiently remote 
from the branch point of the a-expanded eigenval¬ 
ues obtain for A > 0.1. From this we estimate that 
simple analytic expressions for the oscillation proba¬ 
bility accurate to a few percent are quite likely over 
this region using our Hamiltonian formulation 9j. 

To identify the regions where expressions for the 
oscillation probability of our Hamiltonian formula¬ 
tion Q are improvements over simplified oscillation 
probabilities found by Freund. Q, we need to exam¬ 
ine the accuracy of Freund’s full and patched oscil¬ 
lation probabilities. Figures presented in Sect. IVIII 
assess this accuracy. Because the accuracy of Fre¬ 
und’s patched oscillation probabilities is found to be 
vastly inferior to that of his full oscillation proba¬ 
bility, his patched oscillation probabilities are not 
further considered. 

The full a-expanded oscillation probability of Fre¬ 
und Q is compared to the exact oscillation proba¬ 
bility over the interval 0.1 < A < A 2 in Fig. [5] of the 
present paper. We see from this figure that Freund’s 
full oscillation probability [<§] becomes increasingly 
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less accurate as A — > 0 . 1 . On the other hand, the 
oscillation probability of our Hamiltonian formula¬ 
tion evaluated with a-expanded eigenvalues, which 
is compared to the exact oscillation probability over 
the same interval in Fig. 4 of Ref. [9| , remains quite 
accurate as A —y 0.1. We conclude from these results 
that simple expressions for the oscillation probabil¬ 
ity more accurate than Freund’s full result over this 
interval may be found using the results our Hamil¬ 
tonian formulation. 

Freund’s full oscillation probability 0 is com¬ 
pared to the exact oscillation probability in the 
vicinity of the atmospheric resonance in Fig. [ 6 ] of 
the present paper. We see there that Freund’s full 
oscillation probability 0 is discontinuous at the lo¬ 
cation of the atmospheric resonance, whereas the os¬ 
cillation probability obtained from our Hamiltonian 
formulation evaluated with a-expanded eigenvalues 
given in Fig. 6 of our original paper [§] is contin¬ 
uous throughout the atmospheric resonance region. 
Because it is also a very accurate representation of 
the exact result, we conclude that simple expressions 
for the oscillation probability that are more accurate 
than Freund’s full result in the vicinity of the atmo¬ 
spheric resonance may be found using the results our 
Hamiltonian formulation. 

The full a-expanded oscillation probability of Fre¬ 
und 0 is compared to the exact oscillation probabil¬ 
ity for asymptotic values of A in Fig. [ 7 ] of the present 
paper. We see from this figure and Fig. 7 of our 
original paper 0 that Freund’s full oscillation prob¬ 
ability [ 6 | and that obtained from our Hamiltonian 
formulation evaluated with a-expanded eigenvalues 
are of comparable accuracy for asymptotic values of 
A. We conclude that simple expressions for the os¬ 
cillation probability of accuracy comparable to Fre¬ 
und’s full result may be found using the results our 
Hamiltonian formulation for asymptotic values of A. 

In our original publication 0, we also compare ex¬ 
act oscillation probabilities to those found from our 
Hamiltonian formulation evaluated with sin 2 $ 13 - 
expanded eigenvalues. The accuracy with which 
the oscillation probability of our Hamiltonian for¬ 
mulation evaluated with sin 2 (^-expanded eigenval¬ 
ues reproduces the exact oscillation probability is 
strikingly good over the region 0.1 < A < A 2 , as 
seen in Fig. 4 of our original paper. This estab¬ 
lishes that simple and accurate analytic expressions 
for the oscillation probability are also likely to be 
found somewhat above the solar resonance region us¬ 
ing results of our Hamiltonian formulation [9( with 
sin 2 0 i 3 -expanded eigenvalues. 


B. Solar Resonance Region, 0 < A < 0.1 

Similar considerations apply to sin 2 013-expanded 
oscillation probabilities. Close agreement between 
the exact oscillation probability and those of our 
Hamiltonian formulation evaluated with sin 2 $ 13 - 
expanded eigenvalues is seen over the region 0 < 
A < 0.1 seen in Fig. 3 of our original paper 0. This 
suggests that simpler and most accurate analytic ex¬ 
pressions for the oscillation probability within the 
solar resonance region are likely to be found using re¬ 
sults obtained from our Hamiltonian formulation Q. 

The figures in our original paper 0 also show that 
for values of A closer to the branch point of the 
sin 2 0 i 3 -expanded eigenvalues, the accuracy of the 
oscillation probabilities deteriorate rather rapidly. 

AS before, to identify the regions where expres¬ 
sions for the oscillation probability obtained from 
our Hamiltonian formulation [9| are improvements 
over simplified oscillation probabilities found by 
AHLO [3|, we need to examine the accuracy of the 
full AHLO oscillation probabilities. Figures pre¬ 
sented in Sect. I VIll assess this accuracy. 

We see from the results appearing in Figs. [3] to 
Fig. □ of the present paper that the full sin 2 $ 13 - 
expanded oscillation probability of AHLO agrees 
rather poorly with the exact result over all regions 
except for the far solar resonance region. Within the 
far solar resonance region, the full sin 2 (^-expanded 
oscillation probability of our Hamiltonian formula¬ 
tion 0 is reasonably accurate. 

From these results, we conclude that simple ex¬ 
pressions of accuracy greater than a few percent are 
likely to be found from the results our Hamiltonian 
formulation evaluated with sin 2 ^-expanded eigen¬ 
values. Such oscillation probabilities would then be 
vastly more accurate than those of AHLO |7f. 

IX. SUMMARY AND CONCLUSIONS 

In this paper, we assessed the theoretical errors of 
simplified neutrino oscillation probabilities obtained 
from various approximate formulations. The errors 
of the simplified oscillation probabilities obtained in 
Refs. [ 6 h 8 | are associated, in large part, with branch 
points |9|] of the medium-modified neutrino eigenval¬ 
ues. Currently, there is essentially no understanding 
of the errors of these simplified oscillation probabili¬ 
ties. We also assessed the likelihood that simple an¬ 
alytical approximations to three-flavor neutrino os¬ 
cillation probabilities with accuracy of a few percent 
would be found using our recent Hamiltonian formu- 
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lation Q. We considered flavor-changing transitions 
in the sector within the Standard Neutrino 

Model, the only cases considered in Refs. 

We compared the approximate results over various 
regions of A to identify which gave the best descrip¬ 
tion of the exact results. From this comparison, we 
found that the full results of Freund are satisfactory 
everywhere except for the solar- and atmospheric- 
resonance regions. Those of AHLO describe the 
solar-resonance region well, but they are compli¬ 
cated and rapidly deteriorate outside this region. 
The simple neutrino oscillation probabilities of the 
approximate formulation we propose in Sect. ELLO 
are accurate to a few percent errors within essen¬ 
tially all regions of interest. 

We thus conclude from Sect. IVlfll that our Hamil¬ 
tonian formulation offers the opportunity to obtain 
approximate oscillation probabilities vastly more ac¬ 
curate than those presently available. This result 
is significant because it would obviate the need for 
exact computer simulations under many circum¬ 
stances, thus simplifying analysis and prediction 
at current and future and next-generation neutrino 
neutrino facilities. 

Making use of the guidance provided by the results 
of this paper, in future work [lb] we derive simple 
approximations with few percent errors within es¬ 
sentially all regions of interest using our Hamiltonian 
formulation. 


Appendix A: Oscillation Probability in Our 
Hamiltonian Formulation 


In this appendix, we give alternative expressions 
for the sin 2 6 *i 3 -expanded eigenvalues equivalent to 
those given in Ref. j9]. Separate expressions are 
given for A above and below A = a in terms 
of the ratio of the small parameters of the SNM, 
R p = sin 2 6 * 13/0 ~ 0.7 


1. Alternative Representation for the 
sin 2 #i 3 -Expanded Eigenvalues 

±L 

We introduce the notation E ee to represent the 

2 ~ 1 ~ G 

sin 0 i 3 -expanded eigenvalues Eg for A < a, and E et 

to represent them for A > a. The representation of 

- . . - L - G 
Eg given in terms of Egg and Egg is equivalent to 

the representation Eg given in Ref. [9|. 


a. A < a 


- L 

The following expressions for Egg are found by 
replacing A —> aR s and Cg —> olCt in Eg given in 
Ref. @, 


E ei — 2 "(-^s + 1 ~ Ct) 

+ —— 4:Aq + Aq — 2a + 2A^a) 
ZUt 

— L QJ 

E 92 = + 1 + Ct ) 

~ 2{l-y)C T {2 ~ A °~ a + aCT) 

— 2 ^ ~ 4Aq + Aq — 2a + 2Aqo) 


(Al) 


b. A > a 


- G 

The following expressions for Egg are found by 
replacing A -A a/R s and Cg —> ACt in Eg given in 
Ref. H|, 


C = ^(R. + 1 -C T ) 

A/P a 

+ 2R, {1 A)C T {2R -- A ° R -- aR -- aCT) 

+ V oA P ( 2 - 4Ao + Al - 2a + 2A 0 a) 

E 82 = + 1 + @t) 

A/P A 

- 2R, {1 ~ y)cC R -- A ‘ R -- aR - + aCT) 
- y AA(2-4A„ + A‘-2 a + 2A, a ) 


~G 

E 83 


1 + 


a 2 R p 

RsC-y) ' 


(A2) 


c. Discussion 

Because Eqs. (IA1IA2I) are equivalent representa¬ 
tions of the eigenvalues Eg and their differences, they 
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are accurate to a few percent for A < 0.6. The fol¬ 
lowing points are also worthy of noting. 

The same quantity 

C T = V 1 + R 2 s~ 2 R s COS 0\2 

» 1 , (A3) 

appears in the expressions for the eigenvalues for 
both A < a and A > a. 

~ G 

Above the branch point at y = 1, E ge and E ge are 
related; 

(1) Eg 2 = Eqi\c t —>—Ct and 

(2) Eg 3 and Eg 2 exchange roles. 

Appendix B: Exact Analytic Results for the 
Coefficients w < : rnn ' > [f] 

Full expressions for the coefficients uy"''”* appear¬ 
ing in the expressions for the partial oscillation prob¬ 
abilities of Ref. [§] are given in this appendix. They 
are expressed in terms parameters of the SNM and 
the following functions of them, 

= cos 2 0 i 2 cos 2 023 ± sin 2 023 sin 2 0 12 sin 2 0 i 3 
C , 2 ±) = cos 2 0i2 ± sin 2 0i2 sin 2 0i 3 , (Bl) 

I 


where C 2 ^ was defined earlier in Eq. (??). In the 
SNM, many of the terms in are quite small 

and may be dropped in practice. 

1. Coefficients 

COS* 

A 2 

wfoip-fl = —cos 2 013 sin 2 2012 sin 2 0i 3 (l - a sin 2 0 i2 ) 

3 

( 22 ) a • 2 - 2 /) 
w K 21 =-—sin 2 #i 2 sin t 

COS 2 -, 1 A 10 

2 

W m%-,2 = ~f Sin2 2012 Sin2 0 13 (B2) 

2. Coefficients 


w cos-o = — ^ cos 0i3 sin 2012 sin20i 3 ((l — a ) 2 — A(cos20i3 + a(l — 2cos20i3 sin 2 0i 2 ) 

’ 8 

- 2 a 2 sin 2 9 12 C ( 2 +] )) 

^ cos 013 sin 2012 sin 20 i 3 (1 + a - 2 a 2 sin 2 0 i2 ) - Aw£% 
w rn% = ^ cos 0 13 sin 2012 sin 20 13 (1 - a sin 2 0 i2 ) , (B 3 ) 


w cos-0 = 77 cos 013 cos 2023 sin 2012 sin 20 i 3 ((l - a ) 2 - A(cos 20 i 3 - a{\ - 4 sin 2 0 i2 sin 2 0 i3 ) 

’ O 

+ 2 a 2 sin 2 0 12 C { 2 ~ ] )) 

w cos!i = ^ cos 2023 sin 2012 sin 013 (cos 2 0i3 + a cos 2 0i3 - 2 a 2 C 2 ^) 

w ¥o% = - ^ cos 2023 sin 2012 sin 0 13 (cos 2 0 i3 - aC^ _) ) , (B 4 ) 
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w cos!o = 77 sin O13 sin 20 12 cos 2 6*13 ((1 - a) 2 sin 2 0 23 - ^ 4 (cos 20 i 3 sin 2 0 23 

+ a(cos 2 0i2 — sin 2 0 23 (1 — 4 sin 2 0 32 sin 2 0 23 )) — 2 a 2 sin 2 0 i 2 C'} + ■*)) 
uP 2 s l = ^ sin 20 12 sin 0 13 (sin 2 0 23 cos 2 0 13 + a cos 2 0 13 sin 2 0 23 - a 2 (l - 2C'^ +) ) 

+ A cos 2 0i 3 (l — a sin 2 0i 2 )) 

w cos;2 = _ ^ sin 2^12 sin 01 3 (2 cos 2 013 sin 2 0 23 - a(l - 2 C'^ +) )) . (B 5 ) 

/qo^ • ( 22 ) 

The coefficients Wc OS are obtained from w y CO s by making the replacement sin 0 23 o cos 0 23 and by flipping 
the overall sign. 


3. Coefficients 

w£ 2) = ^(1 — a) 2 sin 2 0 12 sin 2 20 i3 sin 2 0 23 + -Acos 2 0i 3 (sin 2 20 i 3 sin 2 0 23 — 4asin 2 0i 2 sin 2 0i 3 sin 2 0 23 

x (2 — 3 sin 2 0 13 ) — 4a 2 sin 2 0 i2 sin 2 0i 3 (C^ + sin 2 0 23 cos 20 i2 — sin 2 0 23 C^ + sin 2 9 23 C^) 

+ 4a 3 sin 2 9 12 C[ +) C^ +) ) 

Wq.^ = — cos 2 0i 3 (sin 2 0 i 3 sin 2 0 23 — a sin 2 0 12 sin 2 0i 3 sin 2 0 23 — a 2 sin 2 0 i2 sin 2 0i 3 sin 2 0 23 
+ a 3 sin 2 9\ 2 C^) — Aw 

u;g 2 2) = c °s 2 013(sin 2 0 13 sin 2 0 23 — 2a sin 2 0 i2 sin 2 0 i3 sin 2 0 23 + a 2 sin 2 0i 2 c} + g , (B6) 

(13) • (12) 

The coefficients Wq ; are obtained from Wq ’ by making the replacement sin 0 23 o cos 0 23 . 
icgg = ^(1 - a ) 2 sin 2 20 23 cos 2 0i 3 C^ _) 

— cos 2 0i 3 (cos 2 0i 3 sin 2 20 23 sin 2 0 33 — a sin 2 9 13 sin 2 20 23 (1 + sin 2 0 32 — 3 sin 2 0 32 sin 2 0i 3 ) 

— a 2 (sin 2 20i 2 sin 2 0i 3 + sin 2 0 32 sin 2 20 23 (1 — 3 sin 2 0i3)C' 2 _ g 
+ a 3 sin 2 0i 2 (sin 2 20 i2 sin 2 0 i3 + sin 2 20 23 + C 2 ' >2 )) 

= - -(cos 4 0 i 3 sin 2 20 23 — a cos 2 0i 3 sin 2 20 23 C 

— a 2 cos 2 0 i 3 sin 2 20 2 3 C 2 _ ^ + a 3 (sin 2 20 32 sin 2 0 33 + sin 2 29 23 C^ 2 )) 

Wq 2 .^ = i(cos 4 0 i 3 sin 2 20 23 — 2 aC^ cos 2 0i 3 sin 2 20 23 + a 2 (sin 2 20 i 2 sin 2 0i 3 

+ sin 2 20 23 C , 2 _)2 )) , (B7) 


4 ;o } = — -j-(l — a) 2 sin 2 0 32 sin 2 20 33 — -Acos 2 0 i 3 (sin 2 20 33 — 2 asin 2 0 32 sin 2 0 i 3 (l + 3 cos 20 i 3 ) 

— 4 a 2 sin 2 0 32 sin 2 0 i 3 (l + cos 2 0 32 — sin 2 0 32 (2 — 3 sin 2 0 i3 )) + 4 a 3 sin 2 0 i 2 C2 + ^ 2 ) 

Wq}P = cos 2 0i3 (sin 2 0i3 — a sin 2 0 32 sin 2 9 13 — a 2 sin 2 9 12 sin 2 0 i 3 + a 3 sin 2 9 i 2 C t ^ r ' > ) 

- Av>$ 

= — cos 2 0i3(sin 2 0i3 — 2 a sin 2 0 32 sin 2 0 33 + a 2 sin 2 0 \ 2 C^) (B8) 
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wffi = —a(l — a) 1 2 cos 2 0i 3 sin 2 0 23 C[ + ^ — A cos 2 0 i 3 (cos 2 d 33 sin 2 0 33 sin 4 5 6 7 8 9 10 11 12 0 23 

— a sin 2 0 i 3 sin 2 0 23 (1 — 3C^ + cos20i 2 cos 20 2 3 ) 

+ a 2 sin 2 0 i2 sin 2 0 23 (1 — 3 sin 2 0i3)Ci + ^ + a 3 sin 2 0i 2 Ci + ^ 2 ) 

4 2 1 2) = cos 2 0 i 3 cos 2 0 2 3(1 — cos 2 0 i 3 sin 2 0 2 3 ) — a cos 2 0 i 3 sin 2 0 2 3 c{ +) — a 2 cos 2 0 i 3 sin 2 d 23 c[ + ^ + oAC^ 
x (1 — Ci + ^) + A cos 2 0 i 3 (sin 2 0i 3 sin 2 0 23 — 2a sin 2 0 12 sin 2 0 i 3 sin 2 d 23 + a 2 sin 2 0i 2 C'i + ^) 

= — cos2 #13 sin 2 # 23(1 — cos 2 0 i 3 sin 2 0 2 3 ) + 2a sin 2 0 2 3 (cos 2 0i 2 cos 2 0 2 3 — sin 2 0 23 sin 2 d 32 sin 4 0i 3 

- sin 2 0i3 (cos 2 0i2 - sin 2 0 23 )) — a 2 c[ + \ 1 - C^ +) ) (B9) 

(33) • (22) 

The coefficients w g ' are obtained from Wq ' by making the replacement sin 0 2 3 o cos 0 2 3 . 
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